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Abstract 

Existing optimal control protocols for mitigating the effects of relaxation and/or RF inhomogeneity on broad- 
band pulse performance are extended to the more difficult problem of designing robust, refocused, frequency 
selective excitation pulses. For the demanding case of T\ and T2 equal to the pulse length, anticipated signal 
losses can be significantly reduced while achieving nearly ideal frequency selectivity. Improvements in perfor- 
mance are the result of allowing residual unrefocused magnetization after applying relaxation-compensated 
selective excitation by optimized pulses (RC-SEBOP). We demonstrate simple pulse sequence elements for 
eliminating this unwanted residual signal. 
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1. Introduction 

Frequency-selective pulses have widespread util- 
ity in magnetic resonance and have motivated sig- 
nificant efforts towards their design 
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3fl 1371 |38|, [391 |40|, |41|, |42j, |43j, |44j, |4£ 
many useful cases, the resulting methodologies can 
achieve the best approximation to the ideal rect- 
angular profile for spin response as a function of 
frequency offset. 

However, in all of these approaches to pulse de- 
sign, performance criteria that can be included in 
the design protocol are restricted either by analyt- 
ical procedures of highly specific scope or by nu- 
merical methods that are limited by the efficiency 
of the optimizations employed. As a result, pulse 
response is typically optimized only for the nom- 
inal ideal RF pulse values. In addition, although 
the length of pulses required for narrowband ap- 
plications can significantly reduce their effective- 
ness when relaxation times are comparable to the 
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pulse length 3, 4|| , the solution to the problem — 
selective pulses which are less sensitive to relaxation 
effects — can also be demanding for standard design 
methods 0, [ED, [H [H [H [Hj . 

To make these design challenges tractable, the 
space of possible pulse shapes is often reduced by 
forcing the solution to be a member of a particu- 
lar family of functional forms (for example, finite 
Fourier series). Thus, potentially, there are impor- 
tant solutions that are missed. 

Over the past decade, we have shown opti- 
mal control theory to be an efficient and power- 
ful method that can be applied to a wide range 
of challenging NMR pulse design problems without 
restricting the space of possible solutions 55, 56, 

69, 
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62, 63, 64, 65J, |£& [67 
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band pulses [6 
are simultaneously tolerant to RF inhomogeneity 
and relaxation effects, which we develop further in 
the present work. 



It is capable of desi gnin g broad- 
and selective pulses [ZJ, [75j that 



2. Selective pulse design 

Optimal control (including similar, related op- 
timization procedures) was originally introduced 
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into magnetic resonance for the design of band 
selective pulses, primarily for imaging 76j, [77 
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[. It was quickly supplanted 
by the efficient Shinnar-LeRoux (SLR) algorithm 
H, IH, S3, S3) which establishes a correspon- 



dence between frequency-selective pulse design and 
digital filter design. There are fast, non-iterative 
algorithms for the ideal filter and, hence, the ideal 
pulse. Unfortunately, the algorithm docs not ac- 
commodate additional criteria, such as tolerance to 
RF inhomogeneity (included in some of the earli- 

~|) or 
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est optimal control-related approaches 
relaxation effects. In addition, the most applica- 
ble and widely used form of the algorithm derives 
pulses which produce a specific linear phase dis- 
persion in the spectral response. Pulses producing 
no phase dispersion, suitable for spectroscopy, are 
more problematic for the SLR algorithm. 

We first provide an overview of well-known issues 
relevant to selective pulse design, since there is con- 
siderably less freedom in the choice of parameters 
compared to broadband pulses. For example, in de- 
signing broadband pulses, we have shown [60] there 
is at best only a marginal relation between the max- 
imum amplitude, -R-Fmax, of a pulse and the achiev- 
able excitation or inversion bandwidth, as long as 
the pulse length, T p , is allowed to increase suffi- 
ciently. Alternatively, increasing RF max for a given 
T p can improve performance for a given bandwidth 
or increase the bandwidth. There can also be in- 
numerable broadband pulses that provide approx- 
imately equal performance for a given RF max , T p , 
and bandwidth. 

Selective pulses are far more constrained, with 
a well-known relation between the selective band- 
width and T p , and much tighter limits on the choice 
of i?T max for a given product of bandwidth and T p 



21| . We provide only a simple overview of the op- 
timal control methodology, emphasizing the modi- 
fications necessary for the present work. The basic 
algorithm for optimizing pulse performance over a 
range of resonance offsets and RF inhomogeneity is 
described fully in 1571 . Details related to incorpo- 
rating relaxation [66| and specific dispersion in the 
phase of the final magnetization [6 ll I67L |83[, which 
we refer to as the phase slope, are provided in the 
associated references. 

2.1. Phase slope 

Values of the phase slope, R, at each offset 67 1 
characterize the net phase dispersion that accumu- 
lates during a pulse of length T p . The phase slope 



is defined relative to the net precession of trans- 
verse magnetization that would be produced solely 
by chemical-shift evolution during the same time 
interval, T p . A pulse that produces focused mag- 
netization of fixed phase for all spins in the offset 
range of interest would have constant R — (i.e., 
a self-refocused pulse). Many selective pulses are 
symmetric, R = 1/2 pulses 0, SI SH, SI SI Sfl| • 
The symmetry of the resulting pulse provides an 
advantage in the development of various algorithms 
used in selective pulse design, such as SLR, inverse 
scattering 
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27J , polychromatic [29|, and stereo- 
graphic projection [30|. In fact, the standard form 
of the SLR algorithm [13 can only generate linear 
phase of this value. Sophisticated algorithms allow- 
ing for more general phase in selective pulses are 
described in the literature 2(| 84 , 8j| , but they are 
specific to this particular performance factor and 
cannot include tolerance to variations in other ex- 
perimentally important parameters. 

By contrast, including additional performance 
criteria, such as a general phase response, is 
straightforward for optimal control. Initial mag- 
netization M(to) is driven by the RF controls to a 
final magnetization F that is defined for each res- 
onance offset in the desired range. To excite trans- 
verse magnetization of linear phase slope R, we con- 
sider target states for each offset w in the excitation 
bandwidth of the form [671 



F = [cos(ip), sm(ip), 0] 



(1) 



Choosing if — RloT p gives a linear phase slope, but 
any function can be considered to define a useful 
target phase, such as quadratic or higher order. 

Selective excitation most simply requires chang- 
ing the target to F = [0,0,1] for offsets outside 
the desired bandwidth. In principle, this stopband 
includes an infinite range of frequencies that must 
therefore be truncated at some chosen value. We 
found as a practical matter that choosing the stop- 
band to be ~ 5 times the passband width was suf- 
ficient to eliminate excitation at higher frequencies 
for the pulse parameters used here. This value can 
easily be adjusted upwards if necessary, or, alter- 
natively, high-frequency components of the result- 
ing pulse determined from Fourier analysis can be 
deleted after verifying they have no significant ef- 
fect on the passband excitation. 

In addition, since the ideal rectangular offset re- 
sponse can not be excited by a finite-length pulse, 
there must be a transition connecting the excita- 
tion frequencies to the nulled frequencies. This is 
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typically defined in terms of design parameters for 
finite impulse response (FIR) filters. An overview 
of the issues relevant to our optimal control imple- 
mentation follows. 

2.2. Selective pulses as digital filters 

For design conditions employing ideal RF in the 
absence of relaxation, selective pulse performance 
is completely determined by the desired passband 
width B, pulse length T p , transition width W join- 
ing the passband and stopband, and residual signal 
fluctuation or ripple 5i and 82 about the ideal target 
amplitude for the passband and stopband, respec- 
tively. 

The passband frequency v v and stopband fre- 
quency v s are defined where the magnitude of the 
magnetization response becomes less than the as- 
sociated fluctuations 1 — 8\ and j^l, as illustrated 
in Fig.l (adopted from Ref. The frequency 

where the amplitude drops to one-half is approxi- 
mately the average of these two frequencies. The 
full width of the filter is defined as twice this value, 
giving a bandwidth B = v s + v v and a fractional 
transition width W = (y a — v p )/B. 

More specifically (and again emphasizing the de- 
sign conditions stated at the beginning of the sec- 
tion) , selective pulse performance is constrained by 
relations for optimal FIR filters of the form 

WT p B = f(6 1 ,6 2 ), (2) 

in terms of an empirically derived function f(5%, 82) 
[jij . For a given value of / = WT p B, smaller 
(larger) 81 gives larger (smaller) 62- Alternatively, 
for fixed 5\ or fixed 82, values of / increase as 62 
or S\, respectively, decrease. Flexibility in selec- 
tive pulse design is thus purchased at the cost of 
trade-offs among bandwidth, pulse length, transi- 
tion width, and ripple amplitudes. Choosing any 
four of the set determines the fifth. 

This relation appears to have been little used in 
the spectroscopic community. Although the pre- 
cise form of the function f {81,82) holds only for 
R = 1/2 pulses, we have found it to be a useful 
qualitative indicator for more general R. One im- 
portant implication is that pulse performance for a 
given absolute transition width BW = v s — v p can 
be made independent of the passband width, B. 
Fixed T p results in the same performance in terms 
of residual signal (ripple) for different B as long as 
the transition width BW is constant. This was ob- 
served empirically and noted in [47[. We thus use 
Eq. to inform our optimal control design. 



2.3. Optimal control 

The approximation to the ideal rectangular fre- 
quency response profile, as illustrated in Fig.l, can 
be readily obtained using a variety of methods — 
among them, the Parks-McClellan (PM) algorithm 
for linear-phase FIR digital filters [S7J ■ After calcu- 
lating the PM polynomial, the standard SLR algo- 
rithm effectively inverts the frequency response to 
produce a linear phase (R = 1/2) RF pulse. The 
equivalent optimal control approach would be to 
use the polynomial response function as the target 
response and efficiently modify the RF controls to 
achieve the allowed performance. 

However, as noted already, the target response 
derived using the PM algorithm applies only to 
R = 1/2, using the ideal RF amplitudes, in the ab- 
sence of relaxation. More significantly, optimal con- 
trol does not need to know the polynomial response. 
One can simply define the ideal passband/stopband 
frequency response, and the algorithm will find the 
response allowed by the particular choice of band- 
width, pulse length, and transition width. Different 
functional forms can be defined for the response in 
the transition region to provide additional flexibil- 
ity. The response at each frequency can also be 
given weights to fine-tune the final excitation pro- 
file. 

We now proceed with the design of more robust 
selective excitation pulses. In what follows, we let 
the frequency response to the pulse (ie, the tar- 
get function) be undefined in the transition region. 
The pulse has the flexibility to do anything it wants 
there. An important addition to the algorithm is 
an adjustable weight function which changes the 
weight for each particular offset depending on the 
deviation of intermediate results from the desired 
performance. During a given iteration, if the de- 
viation of magnetization from the target for a par- 
ticular offset is larger than the allowed ripple am- 
plitude, then the weight for this offset increases; 
otherwise it decreases. This is an extension and 
generalization of the method introduced in [6l| . 

3. Results and Discussion 

We present several examples illustrating the ad- 
vantages of relaxation-compensated selective exci- 
tation by optimized pulses (RC-SEBOP). As dis- 
cussed, possible performance improvements include 
increased tolerance to RF inhomogeneity and com- 
pensation for relaxation effects. 
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3.1. Compensation for relaxation and RF inhomo- 
geneity 



The algorithm for generating relaxation compen- 
sated broadband pulses (RC-BEBOP) has already 
been developed [66|. It only needs inclusion of 
the modified target functions described in Section 
I2.1l for application to selective pulses (RC-SEBOP). 
The primary result of that work is that substan- 
tial signal gains are possible relative to expected 
losses from short Ti, T 2 by finding trajectories for 
the transformed magnetization that minimize these 
losses, even if this requires a longer pulse length. 
Relaxation losses are minimized by keeping spins 
near the z-axis and orienting them so that all off- 
sets can be transformed to the transverse plane in 
a very short time at the end of the pulse. This so- 
lution not only mitigates the effects of transverse 
relaxation, but short Ti becomes an advantage due 
to repolarization of z-magnetization during the rel- 
atively long time the spins are aligned close to the 
longitudinal axis. The net affect is that almost all 
the RF power is applied at the end of our relax- 
ation compensated pulses [66[. This is fortunately 
also consistent with a particular family of solutions 
for broadband R = (refocused) pulses 67 1. 

There are also solutions for refocused selective 
excitation pulses [i?} that employ significant RF 
power throughout the pulse and therefore do not 
lend themselves to relaxation compensation. On 
the other hand, these pulses do an excellent job of 
minimizing the residual off-axis magnetization in 
both the passband and stopband. Consistent with 
this result, we find empirically that minimizing re- 
laxation losses for selective pulses competes with 
minimizing residual off-axis magnetization. The 
trajectories that reduce relaxation effects are not 
compatible with those that achieve good refocus- 
ing. We therefore consider a strategy that maxi- 
mizes a;-magnetization in the passband, minimizes 
it in the stopband, i.e., the standard procedure, 
but removes any explicit restrictions on residual y- 
magnetization. This allows the optimal control al- 
gorithm to emphasize relaxation compensation and 
find solutions with considerably enhanced perfor- 
mance compared to those that give highest priority 
to minimizing the y-componcnt. As we will show, 
this undesirable residual magnetization can be elim- 
inated after the pulse without significantly affecting 
performance. 



3.1.1. Excitation profile effects 

In addition to causing signal losses, relaxation 
can also degrade the uniformity of the excitation 
profile [II]. SLURP pulses were developed 
specifically to obtain more uniform response over 
the selected bandwidth while accepting attenuation 
due to short Ti, T 2 . SLURP- 1 pulses address the 
particular case T\ = T 2 , and were derived for vari- 
ous values of the ratio T 2 /T p . 

For the demanding case T p = Ti = T2 = 1 ms, 
theoretical performance of SLURP- 1 is compared 
to RC-SEBOP in Fig. for RF miscalibrations 
of -10%, 0%, and +10%, displayed left to right. 
For the ideal RF shown in the middle panel, RC- 
SEBOP provides a signal gain of 60%, uniformly 
preserving greater than 90% of the x-magnetization 
over the optimized bandwidth of 4 kHz, in spite of 
short Ti,T 2 equal to the pulse length. It is rela- 
tively insensitive to RF inhomogeneity of ±10%. In 
all cases, there is improved signal suppression from 
M x magnetization in the stopband, and a narrower 
transition width. 

Analogous investigations have optimized 
relaxation-compensated pulses by applying simu- 
lated annealing to RF waveforms represented using 
finite Fourier series (~ 8 frequency components or 
less) and no compensation for RF inhomogeneity 
[HI HI. They achieve signals of < 80% for the 
case T p = T\ = T 2 over a less selective bandwidth 
than obtained here. Some of our performance 
gains may be due to a more efficient optimization 
procedure that does not restrict the solution space. 
However, we expect the largest gains are due to 
the flexibility of allowing residual y-magnetization, 
My, which can be quite large, as shown in the 
bottom panel of Fig.j^l This is easily eliminated, as 
shown later, and allows a much more ideal, rectan- 
gular excitation profile than previously considered 
possible for short Ti, T 2 481. with minimal loss of 
signal. Alternatively, if we choose to minimize the 
residual M y during the optimization, this requires 
trajectories that sacrifice signal enhancement and 
selectivity of the passband, consistent with the 
other cited studies. 

3.1.2. Excitation using small time-bandwidth prod- 
uct pulses 

A simple strategy for reducing relaxation effects 
is to reduce the pulse length, but according to 
Eq. [[5]] , achieving acceptable performance for a nar- 
row excitation bandwidth then becomes more prob- 
lematic. A given time-bandwidth product requires 
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trade-offs in the sharpness of the excitation profile 
(transition width) and the signal variation (ripple) 
in both the passband and the stopband. Still, one 
can optimize performance for a desired low value 
of the time-bandwidth product, as accomplished in 
the SNOB family of pulses [33|. Including relax- 
ation compensation beyond what is accomplished 
by a short pulse length alone and including toler- 
ance to RF inhomogeneity provide additional op- 
portunities for improved performance. 

For e-SNOB, T p = 1 ms, selective bandwidth 
±1.5 kHz, we designed a RC-SEBOP pulse incorpo- 
rating relaxation times T\ — T2 — T p and tolerance 
to RF inhomogeneity of ±10%. Pulse shapes are 
shown in Fig.[3] along with theoretical performance. 
Signal gains of a factor of 2 are obtained with the 
relaxation-compensated pulse, with a sharper exci- 
tation profile and improved suppression of M x mag- 
netization in the stopband. To minimize relaxation 
losses, RC-SEBOP delivers most of its RF power 
at the end of the pulse. As in the previous exam- 
ple, these enhanced performance features are pur- 
chased at the cost of a larger residual M y compared 
to e-SNOB. In the next section, we present meth- 
ods for selecting only the desired M x magnetization 
while maintaining the performance advantages of 
RC-SEBOP. 

4. Experimental 

All the selective pulses considered so far, both 
traditional pulses and optimal control pulses, pro- 
duce significant residual M y magnetization in the 
passband at non-ideal RF calibration. Methods for 
removing this unwanted signal therefore have more 
general applicability. 

To destroy undesired M y after selective excita- 
tion in a single-acquistion experiment, a hard 90°_ y 
flip-back pulse can be applied to store M x along 
the z-axis. A gradient pulse is then employed to 
dephase the M y component, followed by a hard 
90° pulse for acquisition of the signal due to M x . 
This "crusher" sequence, implemented in Fig. 0^ 
as a more general phase-cycled sequence (to be 
discussed below), was first tested for the case of 
unlimited-bandwidth hard pulses by applying the 
hard pulses on resonance. The transmitter fre- 
quency was shifted only for the RC-SEBOP pulse 
to measure its off-resonance performance. This se- 
quence is insensitive to the actual T2 of the sample, 
since there is minimal relaxation during the short 
hard pulses and no T2 relaxation of magnetization 



stored along the z-axis during the gradient pulse. In 
addition, for typical samples with T\ 3> T2, there 
are minimal T\ effects as a result of the sequence. 
However, for very short Ti, repolarization of stop- 
band magnetization during the gradient pulse can 
lead to slightly more M x magnetization in the stop- 
band than expected from the theoretical selectiv- 
ity profile of RC-BEBOP, previously illustrated in 
Fig.H 

To eliminate extra stopband signal in cases where 
T\ is too short for ideal performance of the original 
crusher sequence, it can be phase-cycled as shown 
in the figure. The two acquisitions add construc- 
tively in the passband. Repolarization at stopband 
frequencies leads to ±z-magnetization after the first 
hard pulse in each cycle which is canceled by addi- 
tion of the two acquisitions after the second hard 
pulse. As a general strategy, this sequence works 
very well for the case of ideal hard pulses with no 
bandwidth limitations. It therefore also works well 
over a bandwidth of ~ 8 kHz for a 25 kHz hard 
pulse, where the performance of the hard pulse is 
sufficiently ideal, and may be useful for applications 
that require only a relatively narrow range of stop- 
band frequencies compared to the available RF am- 
plitude of the hard pulse. 

More realistically, all pulses have to be applied at 
the same transmitter frequency to measure the per- 
formance of the sequence at larger resonance offsets. 
The effective field of the hard pulses will then have 
a significant z-component, and the y- magnetization 
will no longer remain untouched. However, the gra- 
dient still dephases the transverse magnetization re- 
maining after the first hard pulse. To a first approx- 
imation, there is only z-magnetization prior to the 
second hard pulse, and the phase cycle produces 
transverse components of opposite sign that cancel 
on addition. 

Experimental off-resonance performance of RC- 
SEBOP and the phase-cycled crusher sequence is 
demonstrated in Fig. [5] for the residual HDO sig- 
nal in a sample of 99.96% D2O, doped with CUSO4 
to relaxation times of T\ = 1.345 ms and T2 
1 .02 1 ms at 298° K. 

Signals representing offsets between —15 kHz to 
15 kHz were obtained in steps of 200 Hz by offset- 
ting the transmitter in this fashion for all applied 
pulses. The sequence was first implemented using 
10/is hard pulses, which resulted in stopband signal 
of ~ 5% relative to the centerband. It was then 
fine-tuned by increasing the hard-pulse length to 
10.68/xs in the first cycle. The resulting selectivity 



5 



profile is in very good agreement with the simula- 
tions for eSNOB and RC-BEBOP performance in 
Fig-El We obtain a signal of 0.83 on resonance using 
the nominal ideal RF values for RC-SEBOP com- 
pared to the theoretical value 0.89. All values are 
relative to an ideal signal of 1 for the case of no 
relaxation. The small disagreement between exper- 
iment and simulation can be attributed to RF in- 
homogeneity/miscalibration in the flip-back pulse, 
which will leave some small fraction of the desired 
passband signal in the transverse plane to be de- 
stroyed by the gradient. 

Since relaxation compensation functions by keep- 
ing spins close to the z-axis as long as possible dur- 
ing the pulse, the performance of RC-SEBOP is not 
highly specific to particular Ti, T-i values. Although 
optimized for T\ = T2 = 1 ms, the T p = 1 ms 
RC-SEBOP pulse of Fig. [3] performs well for much 
shorter relaxation of T\ — 708 fis and T2 = 527 fj,a, 
as shown in Fig. [6] A signal of 0.73 is obtained on 
resonance using the nominal ideal RF values for the 
pulse and the phase-cycled crusher sequence com- 
pared to a theoretical value of 0.81 for on-resonance 
excitation by RC-SEBOP alone. The pulse could 
also be optimized for faster relaxation to improve 
performance further. 

If a single-acquisition sequence is preferred, the 
sequence of Fig.2}} can be used to more completely 
eliminate stopband signal for the case of short T\. 
The first gradient is followed by a hard 180° pulse. 
M z is flipped to — z where any additional magne- 
tization repolarized during the first gradient pulse 
can relax to zero during a subsequent delay of the 
same length as the first. This extra delay can also 
be used for additional dephasing of unwanted trans- 
verse magnetization by a gradient pulse of opposite 
sign to the first gradient. The sequence ends with 
a hard 90% pulse followed by signal acquisition. 
However, performance depends more sensitively on 
any RF inhomogeneity or miscalibration of the hard 
pulses, since there are now two opportunities to de- 
stroy leftover transverse magnetization due to im- 
perfect rotation and storage along the z-axis. Off- 
resonance effects of the hard pulses, which now in- 
clude a 180° pulse, further degrade performance. 
For the experiment described in Fig. [51 a signal 
of 0.68 is obtained using the single- acquisition se- 
quence of Fig.HjD. The scheme is included as an op- 
tion for narrow passband applications. In addition, 
one might want to explore possibilities using broad- 
band, shaped 90° and 180° pulses. These could also 
be incorporated into a simpler sequence that elimi- 



nates the gradient in FigJUi and cycles according to 
[±y, y, Acq(x, — x)}. These topics are beyond the 
scope of the present article. 

5. Conclusion 

Pulses which provide robust and enhanced per- 
formance despite RF inhomogeneity/miscalibration 
and relaxation effects are highly desirable. The op- 
timal control approach to designing refocused se- 
lective excitation pulses with these compensatory 
mechanisms has been presented. The examples con- 
sidered the standard selectivity profile comprised of 
a passband, transition region, stopband, and vari- 
ations (ripple) in signal uniformity, as illustrated 
in Fig.Q] Constraints and trade-offs in the perfor- 
mance among these fundamental parameters were 
emphasized. In particular, we found relaxation 
compensation and null excitation of magnetization 
in the stopband to be competing goals. Consid- 
erable improvements in selectivity and relaxation- 
compensation for short T\ , T2 were obtained by al- 
lowing residual unrefocused magnetization in both 
the passband and stopband. This residual mag- 
netization can be readily eliminated without sig- 
nificantly diminishing performance using additional 
pulse elements. 
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Figure 1: (Adapted from Ref. [2llH A finite length selective pulse can only approximate the ideal rectangular frequency response. 
Residual signal or ripple amplitude in the selected frequency spectrum (passband) is denoted by Si, with &2 representing the 
ripple over the frequency range where the signal should be nulled (stopband). The positive frequencies v v and v B define the 
passband and stopband, respectively. The plotted response is symmetric about the zero frequency. 



10 




Figure 2: Excited magnetization M x (top panels) and M y (bottom panels) is plotted as a function of resonance offset for 
SLURP-1 (green) and RC-SEBOP (blue) for T p = Ti = T 2 = 1 ms and RF inhomogeneity/miscalibration of -10%, 0%, +10% 
reading left to right. Despite these fast relaxation rates, RC-SEBOP preserves > 90% of the desired M x while achieving a 
nearly rectangular profile that is relatively insensitive to ±10% RF inhomogeneity over the optimized excitation bandwidth 
of 4 kHz. The signal gain on resonance for the nominal ideal RF is a factor of 1.6. Minimal relaxation losses are achieved 
by allowing large residual M y , particularly in the stopband. This unwanted signal can be subsequently eliminated without 
significantly affecting performance, as described in the text and FigsHHBI 
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Figure 3: Selective pulses eSNOB (green) and RC-SEBOP (blue) of length T p = 1 ms arc plotted in panel A. RC-SEBOP 
utilizes less RF power, which is applied primarily at the end of the pulse to reduce relaxation losses, as described in the 
text. Panels B, C, and D show the frequency response of the pulses for RF miscalibaration/inhomogeneity of —10%, 0%, and 
+ 10%, respectively, for T\ = T2 = 1 ms. RC-SEBOP significantly reduces relaxation losses and provides a sharper and more 
rectangular selectivity profile over the designed excitation bandwidth of 3 kHz. The signal gain in this example is a factor of 2. 
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Figure 4: (a) Phase-cycled sequence consisting of excitation by selective pulse (SP), hard 90° pulse, gradient pulse, hard 90° 
pulse, and acquisition. The sequence is designed to eliminate residual ^/-magnetization produced either by standard selective 
pulses in the presence of RF inhomogeneity or produced by design in the case of RC-SEBOP. If Tj is sufficiently long so there 
is no repolarization during the 90°— G — 90° sequence, the first cycle is not needed, and the second cycle can be employed as a 
single-acquisition sequence, (b) a single-acquisition sequence designed as an alternative to (a) for the case of short T\. Further 
details for both sequences are provided in the Experimental secion. 
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T p = 1 ms, T,= 1 .345 ms, T 2 = 1 .024 ms 
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Figure 5: Experimental selectivity profiles of RC-BEBOP (right) and eSNOB (left), pulse length T v = 1 ms, using the phase- 
cycled crusher sequence of Fig. [4^, applied to a strongly relaxing sample with T\ = 1.345 ms and T2 = 1.024 ms. Further 
experimental details are in the text. Results are shown for three values of the RF calibration relative to the ideal value B® 
for each pulse, showing inscnsitivity of RC-BEBOP to RF miscalibration of ±10%. Passband signal gains of a factor of 2 are 
obtained with RC-BEBOP compared to eSNOB. The combined sequence for relaxation compensation and elimination of My 
preserves 83% of the initial magnetization on-resonance for the case Bi/B® = 1 compared to the theoretical value of 89% for 
M x alone shown in Fig. [3] 



14 



T p = 1 ms, T,= 0.708 ms, T 2 = 0.527 ms 
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Figure 6: Same as Fig. [5] but applied to a faster-relaxing sample with T± = 708 fis, T2 = 527/^s. Although RC-SEBOP was 
optimized for T\ = T2 = T p = 1 ms, it provides excellent resistance to relaxation for much shorter values, preserving 75% of 
the initial magnetization on-resonance for the case B\jB\ = 1. 
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